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I. INTRODUCTION 

A combined theoretical analysis of uo and meson production in the processes 
pp — > ppV , pn -^ pnV and pn -^ dV (here V denotes a vector meson [uo or 0], p {n) 
denotes a proton (neutron), and d stands for the deuteron in the final state) at near- 
threshold energies is interesting for different aspects of contemporary particle and nuclear 
physics. For instance, according to the Okubo-Zweig-Iizuka (OZI) rule [1] the production 
of mesons in nucleon-nucleon collisions should be strongly suppressed relative to u pro- 
duction. An enhanced production would imply some exotic (e.g., hidden strangeness) 
components in the nucleon wave function. The OZI rule is based on Sakurai's observation 
[2] that the lowest 1^ vector mesons obey the Gell-Mann SU(3) octet classification [3] and 
the Gell-Mann-Okubo mass formulae only if one attaches to the eight SU(3) matrices Ai_8 



a ninth one, Ag = i/2/3Ao, so that instead of the 1~ octet one considers a nonet, repre- 
sented by a non-traceless tensor G^y ~ q^qi, {qi = u, d, s). Then, to reconcile the physical 
masses of u and mesons one introduces a mixing angle 6 and forms combinations like 
cos^ ujo ± sin^ 0^8 to reproduce the known masses (cjq.s are the pure SU(3) uj meson 
states). Alternatively, one can determine the mixing angle from the demand to reproduce 
the quark content of uj {uu + dd) and {ss) mesons. The angle Oq determined from this 
condition is called the ideal mixing angle which is slightly different from 9 obtained from 
mass formulae. Such a difference means that in principle the meson can contain a small 
portion of non-strange quarks and, vice versa, the wave function of uo can contain some 
hidden strange components. In spite of the fact that the nonet classification does not 
have a strict symmetry nature, it has been found to excellently describe the light vector 
mesons. However, such a "nonet hypothesis" needs to be complemented by the restriction 
that in expressions for physically observed processes the trace G^n will never arise (see 
for details, e.g., ref. [4]). 

At the level of quark diagrams this restriction means that topological diagrams with 
disjoint parts of quark lines ("hairpin" diagrams) must be zero (in the ideal case) or 
highly suppressed. This is known as the Quark-Line Rule or OZI rule. In particular, the 



quantities 



^ozi 



^/2T{A + B ^X+{ss)) 



T{A + B ^X + (uu)) + T{A + B-^X + (dd)) 



Xi) 



and 



0. 



ozi 



V2T{A + B -^X + 



Zozi + tan(9 - Oq) 



1 - Zozi tan 



1.21 



T{A + B ^ X + iu) 

[A, B and X denote non-strange particles, T is the amplitude of the corresponding pro- 
cess) are predicted to be small, 
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The prediction (1.3) has been checked for a large number of experimental cross sections 
and is found to be fulfilled with high accuracy [4], 



f3^ozi < 2.6 ■ 10' 



;i.4) 



Nevertheless, the OZI rule being in fact of a mnemonic nature cannot be exact and should 
have some range of applicability. So, the unitary condition for the S matrix, SS^ = 1, 
implies that 
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where the summation over intermediate states X runs over all open physical channels 
allowed by energy conservation. From (1.5) one concludes that at high enough energies 
OZI allowed subprocesses i — > X and f -^ X always exist such that they can contribute to 
the OZI forbidden i —* f reactions. These correspond to so-called loop or double hairpin 
diagrams being topologically equivalent with the forbidden ones. This means that eq. (1.5) 
leaves room for OZI rule violation. The experimentally observed small values of f3ozi in 
eq. (1.4) may be understood as a random cancellation of intermediate phases in (1.5) at 
high energies, so that the OZI rule can be still fulfilled at high energies [43], while at 
intermediate energies its violation is always expected. Nevertheless, at low energies, near 
the threshold, the loop diagrams contributing in (1.5) can be suppressed by the lacking 
energy, and the hairpin diagrams again govern the amplitude of the process. Hence, since 



in this kinematical region there are no other "legal" sources of the OZI rule violation, an 
investigation of processes near the threshold is of interest. 

Any significant deviation oi /Sozi from (1.4) near the threshold indicates some "exotics" 
(as the mentioned hidden degrees of freedom) in the wave functions of the involved parti- 
cles. Of particular interest is the presently studied uj and production in nucleon-nucleon 
reactions since the violation of the OZI rule could drastically change the interpretation 
of the quark content of nucleons. Nowadays, for the relevant coupling constants one has 
the following predictions (cf. ref. [3]) 

g<t>NN = — ta.n{9 — 9q) Quinn-, 

%^~(0.7-l.)-10-2, (1.6) 

where A0 = ^ — ^q — 3.7° ■■■5.5°. Consequently, within a simplified treatment of the 
reaction mechanism, the ratio of the corresponding cross sections is expected to be pro- 
portional (with the proportionality coefficients corrected by corresponding phase space 
volumes) to the ratios of the coupling constants in eq. (1.6), so that possible violations of 
the OZI rule are often associated with these values [5, 6]. However, the reaction mecha- 
nism of A^A^ -^ NNV, where A^ denotes the nucleon, is much more involved and consists 
of different types of diagrams with quite complicate interference effects. This hinders a 
direct investigation of the validity of the OZI rule; some enhancements of the ratio (1.4) 
may occur dynamically, i.e. the actual ratios of the cross sections may differ from the 
"OZI correct" input ratios of the coupling constants. Moreover, in processes of the type 
pp -^ ppV the effects of Final State Interaction (FSI) may become predominant near the 
threshold and completely mask the studied problem. 

For a reliable study of these effects one needs more experimental data and more types 
of processes. In particular, for further checks of the reaction mechanism and for a firm 
separation of FSI effects it is necessary to study meson production also at neutron targets. 
Near the threshold, FSI in pp and pn systems differs due to the Pauli principle, hence a 
combined study of pp and pn reactions will enlighten the theoretical methods to treat the 
FSI. Unfortunately, data on elementary reactions on neutrons are scarce since they must 
be extracted, with some efforts and even mostly with some model dependent assumptions. 



from reactions on nuclei, mainly on the deuteron. The spectator technique [7, 8] represents 
one example of how one can use a deuteron target to isolate reactions on the neutron. It 
is based on the idea to measure the spectator proton, psp, at fixed beam energy in the 
vector meson production reactions pd -^ dVpsp-, thus exploiting the internal momentum 
spread of the neutron inside the deuteron. In such a way one gets access to quasi-free 
reactions pn — > dV . 

An experimental investigation of the near-threshold (pseudo) scalar and vector meson 
production at the neutron becomes therefore feasible. Indeed, at COSY the ANKE spec- 
trometer set-up can be used in particular for studying the ag, oo and production with 
the internal beam at "neutron target" [8]. This offers the possibility to enlarge the data 
base on hadronic reactions and to address special issues, e.g., for a systematic study of the 
OZI rule violation via u and production in vrA^ and pp reactions (cf. [9] for a reanalysis) 
and in pp annihilations (cf. [10, 11] and further references quoted therein for theoretical 
analyzes) as well. 

OZI rule violations are of interest with respect to possible hints to a significant ss 
admixture in the proton, as supported by the pion-nucleon sigma term [12, 13] and in- 
terpretations of the lepton deep-inelastic scattering [14]. Besides the impact on hadron 
phenomenology the origin of the OZI rule addresses also a link to QCD [15, 16]. Further- 
more, the effective description of particle production in elementary processes is a necessary 
prerequisite to analyze heavy-ion collisions in detail and to pin down in-medium effects. 
In particular, the pn channels deserve a reliable description which is not simply accessible 
from constant isospin factors correcting the cross sections in pp channels. 

Given this motivation, in [17-19] the reaction pn — > dV with V = u^ip has been studied 
in some detail {pn — * dS with S = aQ,!],!]' is considered in [17, 18]). In [17, 18] the cross 
sections and angular distributions are elaborated as a function of the excess energy within 
a two-step model. The same observables are evaluated in [19-21] within the framework 
of a boson exchange model with emphasis on the ratio of cross sections cFpn^d^i/'^pn'^du) 
being of direct relevance for the OZI rule violation. Ref. [22] focuses on the uo and 
meson production in pp reactions. We go beyond [22] by studying also pn reactions and 
by including the deuteron final state. 



In this paper we present a combined analysis of cj and meson production in pp -^ ppV, 
pn -^ pnV and pn -^ dV processes. Within an effective meson-nucleon theory we compute 
the covariant amphtudes for A^A^ — > NNV, and from pp -^ ppV data we fix the free 
parameters as to obtain a good description of the available experimental data. Note 
that in spite of the large number of parameters entering the effective meson-nucleon 
theory a bulk of them is already determined by other independent considerations (the 
one-boson-exchange (OBE) potential, decays into mesons etc.) so that we are left with a 
restricted number of free parameters which can be varied. The main idea of the present 
work is to study whether it is possible to describe in a consistent way the pp -^ ppV, 
pn -^ pnV and pn -^ dV amplitudes by exploiting as input into the calculations such 
effective parameters which, at the elementary level, are in a concord with other data (e.g., 
known meson and nucleon decays) and preserve the OSI rule. Then from calculations of 
the corresponding cross sections we study the possible enhancement of the respective ratios 
and compare with the expected naive OZI rule predictions. For a consistent treatment of 
all mentioned reactions and in order to be able to use directly the covariant amplitudes 
from A^A^ — ^ A^A^V" processes, with the effective parameters already found, we perform 
the analysis of the pn —^ dV processes within the Bethe-Salpeter (BS) formalism [24] 
with the numerical solution [25] of the BS equation obtained within the same effective 
meson-nucleon theory. The use of the BS formalism is not dictated by the necessity of 
taking into account relativistic effects, rather it is inspired by convenience reasons. 

This paper is organized as follow. In Section II the vector meson production in nucleon- 
nucleon interaction is analyzed. The kinematics, notation and explicit expressions for the 
relevant quantities are presented in Section II A, the choice of the effective parameters is 
discussed in Section II B, and in Section II C results of numerical calculations of angular 
distributions and the energy dependence of the total cross sections for uj and meson 
production in pp and pn reactions are presented. Basing on the obtained results the ratio 
of cross sections (f)/uj is analyzed in connection with the OZI rule. A similar structure has 
the Section III, where the meson production in the pn -^ dV process is analyzed within 
the Bethe-Salpeter formalism. In Section III A details of derivation of the corresponding 
formulae within the BS formalism are presented. In Section IIIB results of numerical 
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calculations of the angular distributions, total cross sections and OZI rule are discussed. 
Conclusions and the summary are collected in Section IV. 

II. THE PROCESSES NN -^ NNu AND NN -^ NN(f) 
A. Kinematics and Notation 

Consider the vector meson production in A^A^ collisions of the type 

Ni + N2^N[ + N!, + V. (2.1) 

The invariant five-fold cross section is 

'''■^= ,y, \ A t. E K.U.'fi'r, 1-^, (2.2) 

2^s(s - 4m2) 4 ,,,,2 s'^,s'^Mv ^■ 

where Sj and M.y are the projections of the nucleon and meson spins on the quantization 

axis, and the factor — accounts for n identical particles in the final state. The invariant 
n! 

phase space volume dr-f is defined as 

rfV = {2nY5{p,+p,-p\-p', - q) ^ft...ft...fl.. - (2-3) 



2Ep;(27r)3 2Ep.(27r)3 2Eq(2 
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In eqs. (2.2) and (2.3) the 4-momenta of initial (^1,^2) and final (^1,^2) nucleons and 



vector meson (g) are p = {Ep, p) with Ep = \Jm^ + p^, q = (E^, q) with E^ = Jniy + q^, 
where m and my are the nucleon and meson masses, respectively. The initial energy 
squared of incident nucleons is defined as s = (pi + ^2)^- It is seen from (2.3) that the 
cross section eq. (2.2) is determined by five independent kinematical variables, the actual 
choice of which depends upon the goals of the attacked problem. In the present paper we 
are interested in studying the angular distributions of the produced mesons in the center 
of mass (CM) of initial particles and the energy dependence of the total cross section. 
For this sake it is convenient to choose the kinematics with two invariants, t = {pi — qY 
and S12 = {p'l + ^2)^! and two angles, the solid angle dVL\2 in the CM of the two final 
nucleons and the azimuthal angle (pv of the meson in the CM of initial particles. (This is 
the Chew-Low kinematics.) Then 

d'a = j , \~ ^fi';' \ E \T^-, , I'dtdsudwdni, ^. (2.4) 



The invariant amplitude T / , , is evaluated within a meson-nucleon theory based on 
effective interaction Lagrangians which includes scalar (a), pseudoscalar (tt), and neutral 
{uj,(f)) and charged (p) vector mesons (see e.g. [26, 27]) 

C.NN = 9aNN^^, (2.5) 

C^NN = -^-^^Ni,i^d^{r^^)N, (2.6) 



C-pNN = -9pNN i^N^pTN^/ - -^Na^^rNd"^/) , (2.7) 

jCvnn = -9VNN [ni^N^^ - ^Na^^Nd"^';,^ ; (2.8) 

jCp^v = 9p.v e^.a/3 5^ '^v Tr (9" ^/^ ^„) , (2.9) 



where N and ^ denote the nucleonic and mesonic fields, respectively, Tr ((pp^^) = ^^ <^^ + 
^j" <?^ + <?p" ^^ and bold face letters stand for isovectors. All coupling constants with off- 
mass shell particles are dressed by monopole form factors Fm = (A^ — M^) / (A^ — A;^), 
where A;^ is the 4-momentum of the virtual particle with mass Mm ( k'^ ^ M^). 

These Lagrangians result into two types of Feynman diagrams: (i) the ones which 
describe the meson production from the processes of one-boson exchange (OBE) between 
two nucleons accompanied by the emission of a vector meson from nucleon lines in VNN 
vertices (in what follows we call these diagrams bremsstrahlung type reactions), and (ii) 
production of vector mesons resulting from a conversion of virtual (exchange) tt and p 
mesons into a real vector meson, i.e., from the pnV vertex defined by eq. (2.9)) (these 
diagrams are called internal conversion type diagrams). 

It is convenient to define for the considered processes an effective scattering operator 
0(12; 1'2'y) which can be derived from the Feynman diagrams for the amplitude by 
cutting the external nucleon spinors, as depicted in fig. 1, where the external nucleons 
are represented by solid truncated lines. All the dependencies on meson and nucleon 
propagators as well as on the polarization of the final meson are included into the definition 
of 0(12; 1'2'\/). Accordingly, the invariant amplitude T^ ^^^; ^, can be obtained by merely 
sandwiching this operator between the nucleon spinors. For example, for the pp reaction 
one has 

^itlL'„4 = (K, 4,P2, 4 I 0(12; 1'2V I pi, s,,p2, S2 ) - antisym. (2.10) 



In the case of pn processes, instead of the antisymmetrization in eq. (2.10), denoted 
by "antisym", one should properly take into account the isospin factors (see for details 
ref. [28]). The derivation of the analytical form of the corresponding bremsstrahlung and 
conversion amplitudes (2.10) from the Feynman graphs in fig. 1 from the Lagrangians 
(2.8) is straightforward but rather cumbersome and we do not present it here. Explicitly 
the amplitudes T^ / ^, ^, can be found elsewhere, e.g., in ref. [28]. 



B. Effective Parameters 

In this subsection we discuss the choice of effective parameters for the bremsstrahlung 
and conversion processes. The bremsstrahlung diagrams (fig. la) consist of two parts, 
the pure OBE exchange and the emission of the meson from the nucleon lines. The 
OBE parameters, masses, cut-off's and coupling constants, are adopted to coincide with 
those of the OBE potential by the Bonn group with the pseudo-vector coupling in the 
71 NN vertex [26]. Then in the bremsstrahlung part of the amplitude (2.10) we are left 
with the parameters of the VNN vertices (vertices with wavy lines in fig. la) for the 
real production of u and 0. In principle, the effective parameters for such vertices can 
differ from the corresponding parametrization of the virtual mesons; consequently in the 
present work they are considered as free parameters to be fitted independently. For the 
conversion type vertices pnV (vertices with wavy lines in fig. lb) the coupling constant 
QpTrw has been fixed from a systematic study of pseudoscalar and vector meson radiative 
decays together with the vector meson dominance conjecture (see, e.g. [19, 27, 29]) 
and is found to be (yfp^t^ ~ 12.9 GeV^-^. ( At this point it is worth mentioning that 
often in the literature the definition of the coupling constants Qpj^v differs from ours 
in eq. (2.9) by a mass factor included into Qpj^v as to make it dimensionless. So, in 
refs. [19, 30] this factor reads as JmpVri^ which results in a dimensionless coupling constant 
Qp^^ ~ 10.0.) The corresponding cut-off parameter, also freely adjustable in the present 
consideration, has been adjusted to reproduce the angular distribution of uj mesons from 
COSY-TOF experiments [31] yielding ApTr^; — 1000 MeV . For the production from 
conversion type diagrams the couphng constant gpT,^ is adjusted to data for the ^ ^> pix 



decay. The value r(0 -^ pn) = 0.69 MeV suggests gp-^^f, = —1.1 GeV^^ [28] (cf. also 
ref. [27]). Then for the conversion diagram the cut-off parameter Ap^^^ is found from a 
combined analysis of the contributions of the conversion and bremsstrahlung diagrams to 
the experimental cross section [5, 6]. Together with the nucleonic cut-off A^ = 720 MeV 
we adopt Ap^^ = 1340 MeV, following ref. [28] (set "B"). 

Next we explain the choice of parameters for the bremsstrahlung part. Since the 
nucleon cut-off parameter, An, has already been fixed, we have now to choose only the 
coupling constant QujNN (and k,^) for the emission of the u meson. The corresponding 
parameters for the meson are connected with symmetry relations with the uj meson and 
we do not consider here them as free parameters. The cut-off for uj meson production as 
well as for the meson is chosen as in the OBE potential [26] {A^^n = ^ojAfAf)- In the 
present work the coupling constant QujNN = 9.0 [19] is chosen. Then, having fixed the uj 
parameters, we find QtpNN = — fi'i.cJAfAf tan A^ ^ —0.8 with A^ ^ 5° [4] (cf. also ref. [28]). 
The parameters used in our calculations are listed in Table I. It should be noted that the 
chosen parameters are in the range quoted by eq. (1.6). 

C. Results for NN -^ NNV 

In our calculations we make use of the explicit expressions for the conversion and 
bremsstrahlung diagrams quoted in ref. [28] . The FSI effects have been calculated within 
the Jost function formalism [34] which reproduces the singlet and triplet phase shifts at 
low energies. In fig. 2 we present results of calculations of the angular distribution of uj 
mesons at the excess energy As^'"^ = 173 MeV (As^'"^ = ^/s — 2mp — my) for proton 
(left panel) and neutron (right panel) targets. The experimental data [35] served to fix 
our free parameters for further calculations of the energy dependence of the total cross 
section and for an analysis of pn -^ dV processes. It can be seen that with the set of 
parameters listed in Table I a good description of the data is achieved. It is also found 
that the contribution of bremsstrahlung diagrams (fig. la) is predominant in both pp and 
pn processes (see fig. 4). The difference in magnitudes for pp and pn processes is due to 
the Pauli principle for the former reactions (integration over dVL\2 in eq. (2.2) is performed 
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TABLE I: Model parameters for the process NN -^ NNV (V = uj,(p). 



Vertex Coupling constant 

ifvNN = I^VdVNN ) 



Cutoff 




(MeV) 




^^pTVU 


980 [19] 


^^pn(j> 


1340 [28] 


AtvAT = 


720 


^ujNN '- 


= 1500 



pirV 



VNN 



MNN" 



TT 



a 



QpTTiO ~ 


= 12.9 GeV-^ 


[19, 


33] 


9p-K(f> ~ 


-- -1.1 GeV- 


1 [20, 28] 


QuiNN 


= 9. 






{Ku; = 


0.5) 






g<t>NN 


= -0.8 [28] 






{h = 


0.5, Ae « 5° 


[4]) 




fnNN 


= 1.0 






9aNN 


= 10. 







^<f>NN = ^u}NN = 1500 



A^NN = 1300 
A^NN = 1800 



u 



QpNN = 3.5 



[Up = 6.1) 



QuNN = 15.85 
{n^ = 0.0) 



A 



pNN 



A, 



•NN 



1300 



1500 



"The OBE parameters could, in principle, be different from the ones in real production of on-mass shell 
mesons from NNV vertices. 

only over one hemisphere of the two protons phase volume) and a possible destructive 
interference of diagrams due to antisymmetrization effects and isospin enhancements for 
the latter (see for details ref. [28]). 

The remaining parameters for have been fixed to describe the data for the angular 
distribution of production at As^^"^ = 83 MeV [5, 6]. Having adjusted our parameters 
in this manner we compute the energy dependence of the total cross section. Results 
of calculations together with the available experimental data are presented in figs. 3, 5 
and 6. The dashed lines represent the contribution of the conversion diagrams alone, the 
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dotted lines are results of calculations of both bremsstrahlung and conversion diagrams 
by neglecting the FSI effects. Eventually, the solid lines exhibit the total cross section 
with FSI taken into account [28]. A fairly good description of the data in a large interval 
of the excess energy is achieved. Contrary to uo production, the contribution of the 
bremsstrahlung diagram for (j) mesons is much smaller, due to the small value of the gtf,NN 
coupling. The nontrivial behavior of the ratio <Jpn~>pnV / <^pp~>ppV exhibited in figs. 4 and 7 
evidence that the vector meson production cross sections known by pp reactions can not 
simply be translated into ones for pn reactions. A comparison of the solid and dashed 
lines in figs. 4 and 7 clearly indicates that near the threshold the FSI in pp and pn systems 
is different. Rather, a profound analysis, like the one presented here, is needed to have 
reliable input for simulating heavy-ion and proton-nucleus collisions. This is particularly 
important for the ongoing experiments with HADES [36] . 

Figures 3 and 6 demonstrate that in the extreme threshold limit the main contribution 
comes from the final state interaction effects. As the excess energy increases the role of 
FSI diminishes, even becoming negligible at excess energy As^/^ > 200 MeV . This is a 
disappointing fact, since as mentioned above, the checks of the validity of the OSI rule 
are preferably to be performed at as low energies as possible (where the loop and double 
hairpin diagrams are suppressed), while we find that at these energies the "net" cross 
section relevant for the OZI rule is completely masked by FSI. However, if the FSI effects 
can be firmly separated, the total cross section may still serve as criterion of the validity 
of the OSI rule. This could be estimated by an investigation the OZI (p/oo ratio at equal 
excess energies, where the effects of difference in phase space volumes is minimized, for 
cross sections with and without including FSI corrections. Figure 8 illustrates that the 
dependence of such ratios upon the energy is rather weak and practically is not affected by 
the FSI effects, which let us argue that the adopted Jost method correctly describes the 
FSI effects. The small difference between ratios with and without FSI can be traced back 
to the different interactions in pp {^Sq configuration) and pn {^Si — ^Di configurations) 
systems near the threshold. The results in fig. 8 persuade that in this kinematical region 
the FSI effects, playing an important role in the absolute value of the total cross sections, 
do not essentially mask the study of OZI rule by ratios of total cross sections. Note that 
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the absolute values of the ratios depicted in fig. 8 clearly indicate a violation of the OZI 
rule. But this is an apparent effect since the relative contributions of bremsstrahlung 
and conversion diagrams and, correspondingly, the interference effects are quite different 
for u and production (cf. figs. 3 and 6). Moreover, even at equal excess energies the 
phase space volumes for and uo mesons are quite different. Sometimes in the literature 
one studies the ratio (p/oo at equal initial beam energies [5]. In this case one has two 
effects with opposite contributions: (i) the phase space volume of the meson is much 
smaller in comparison with the u meson (due to different masses, there is a shift of 
~ 240 MeV in As^/^) which will decrease the ratio, and (ii) since at equal beam energies 
the relative momentum of two produced nucleons is smaller for meson production, the 
FSI corrections are expected to be larger (cf. figs. 3 and 6). In fig. 9 the OZI ratio at equal 
beam energies is presented together with data from [5]. It is seen that these ratios are 
by an order of magnitude smaller than that at equal excess energy and more compatible 
with the naive OZI rule predictions. In the right panel of fig. 9 the OZI ratio is depicted 
without any FSI corrections. A comparison with the left panel and with the range of 
the ratio of input constants (the shaded area in fig. 8) gives some evidence about the 
magnitude of corrections from the phase space volumes solely. 

From the performed analysis one can conclude that an investigation of the OZI rule in 
NN -^ NNV processes near the threshold is feasible, provided one can firmly take into 
account the FSI effects and, consequently, properly calculate the phase space corrections. 
The problem of accounting for the FSI between the two nucleons can be solved by con- 
sidering processes with a definite two-nucleon final state, e.g. the process where the FSI 
in a pn system result in a bound state, the deuteron. 

III. THE PROCESSES pn -^ du AND pn -^ d<p 

One way of testing our assertion on the predicted cross sections for pn reactions is to 
implement them in the tagged quasi-free reaction pn -^ dV . 
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A. Formalism 

Let us consider the reaction 

p +n = d + V, (3.1) 

where, as before, V denotes the cu or vector meson and d the deuteron. The invariant 
differential cross section reads 

i =16..(.-4,„=)4,S«E, l«'"'(».*)l^ (3-2) 

where s is the square of the total energy of colliding particles, t is the square of the 
transferred 4-monientum, Si,S2,Aiv and Aid denote the spin projections on a given 
quantization axis, and T stands for the invariant amplitude of the process (3.1). As in 
the previous section the most general form of the amplitude T is presented in the form 

T,^.r-^^(s,t) = (A-M.|G,e;:^Jl,2), (3.3) 

where C,Mv i^ ^^^ polarization vector of the vector meson. The scattering operator G 
represents a vector in the vector space of mesons and a 16 (S> 16 component object in 
the spinor space of nucleons; the deuteron is described as a 16 component BS amplitude 
$(1,2) which is defined as a matrix element of a time-ordered product of two nucleon 
fields iIj{x) as 

$"^(1, 2) = {D\T [iIj"{1)^'^{2)) |0) (3.4) 

and satisfies the BS equation. 

Suppose that in the considered reactions the off-mass shell effects are negligibly small, 
i.e., the on-mass shell matrix elements of the Lagrangian (2.9) between real nucleons 
(reaction (2.1)) and half-off mass shell matrix elements (reaction (3.1)) are the same. 
Then it can be immediately seen that the operator G^ ^* coincides with the corresponding 
operator of the process (2.1), i.e., O = G^ ^*. The invariant amplitude may be cast in the 
form 

T^sl'^^{s,t) = -tj -^^^tS'^\2')d%{l2-l'2\Mv)u%l)u^2\ (3.5) 
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where $^ (1', 2') is the conjugate BS amphtude in the momentum space, p is the relative 
4-momentum of the nucleons in the deuteron, and u{l),u{2) denote the Dirac spinors for 
the incident nucleons. Recall that the operator 0^'^(12; 1'2' ^M.y) is a scattering operator 
describing the vector meson production in the final state. This operator acts in the spinor 
space of protons and neutrons separately; the upper and lower spinor indices refer to 
protons and neutrons, respectively. The first indices, h and a, form an outer product of 
two columns, whereas the second ones, c and /5, form an outer product of two rows. To 
explicitly establish a relation of the amplitude (3.5) with the corresponding amplitude 
(2.10) of the A^A^ -^ NNV process we find the spinor structure by a decomposition of the 
operator O in each of its indices over the corresponding complete set of Dirac spinors, i.e, 

d^^^(12; l'2\Mv) = y^ E <rp,p.(12; l'2')t.^(lO^;:(l)«?(2)«^,(2'), (3.6) 

*> ' r,r',p,p'=l 

where the coefficient A^,^y ,{12; 1'2') is determined by the completeness and orthogonality 
of the Dirac spinors, Ur{p)ur'{p) = er2m6rr', yielding 

A^,X^^^,il2; 1'2') = erer'e,e,,uUl')u''A^')0%il2; l'2',A^v)<(l)<(2), (3.7) 

where e = +1 for r = 1,2 and e = —1 for r = 3,4. Since the defined operator (3.6) 
generally acts in the nucleon spinor space its matrix elements describe processes with 
anti-particles as well. This is explicitly refiected in the dependence of the coefficients 
A^J,{12] 1'2') on spinors with anti-particle quantum numbers. In our case this depen- 
dence can occur in the process (3.1) only as virtual creation and annihilation of NN 
pairs, which within the BS formalism are allowed through the presence of negative-energy 
P components in the BS amplitude. However, in the considered kinematics the contri- 
bution of P waves is by far smaller than the main S^'^ and D^^ components [37] and 
consequently, in what follows we disregard all the partial BS amplitudes with at least one 
negative energy index, leaving only the S^^ and D^^ waves. Then substituting (3.6) into 
(3.5) one obtains 

T^sl^'^is^t) = j£yj 7^^^.(1''2') E <-,.,(12;l'2')<,(2')^^(lO 



rr'=l 



E / 7Sl«'(2')f 7r'(7f""^^'!(l',2')) A^-.^„,{12;1'2')4{1') 



{2m)\^^J (2 
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= (^ E / (0i<.;.w(12; 1'2')vA'^')^mA^\ 2')t..(l'), (3.8) 

where 7c is the charge conjugation matrix, Vr{2') = uj{2')'~fc, and the new BS amphtude 
'l'_A4^(l',2') = 7c$;k^(1', 2') is written now in the form of a 4 4 matrix and represents 
the solution of the BS equation written also in matrix form. Note that in eq. (3.8) the 
spinor Vr'{2') does not describe an anti-particle, it is rather a consequence of the efforts 
made to cast the BS amplitude in a matrix form (see for details, e.g., refs. [24]) and still 
refers to nucleons. 

As in the previous case of A^A^ reactions, the effective interactions in eq. (2.9) 
again result into two types of diagrams, conversion currents and bremsstrahlung emis- 
sion, which are depicted in fig. 10. Consequently, having already computed the op- 
erator O, it is straightforward to obtain the coefficients A^y.pp,{12] 1'2') in eq. (3.7) 
within the OBE approximation. Note that in case when all particles are on mass shell, 
A^J,.pp,{12]V2') exactly coincides with the amplitude (2.10) of the elementary process 
(2.1), i.e., Af/,.pp,{l2; 1'2') = T;^,!;_^,(12; 1'2'\/). However, in general this amplitude cor- 
responds to a virtual process of vector meson production with two off-mass shell nucleons 
in the final state. 

Since our numerical solution for the BS equation has been obtained in the deuteron 
center of mass [25], all further calculations will be performed in this system. First we 
introduce the relevant kinematical variables as follows: pi and p2 are the four momenta 
of incoming nucleons, p'^ and p'2 stand for the momenta of the internal (off-mass shell) 
nucleons in the deuteron with p = {p[ ~p'^j2^ ^j^ denotes the polarization 4-vectors of 
the deuteron and vector mesons. In this notation the BS amplitudes in the deuteron rest 
system are of the form [37] 

^X^,^(Pi,P2) =Ar(^i + m)i±^e^,(^2-m)05(Po,|p|), (3.9) 

'^^mI^A^V'^ = ^7f (^1 + "^)^^ [^■^'i + 2k^*'^^ " ^2)(Km) j ik - rn)(l)D{po, |p|), 
where ki^2 are on-mass shell 4-vectors related to the off-mass shell vectors p[ 2 as follows 



h = {Ep, p), k2 = {Ep, -p), p[ = {p[q, p), P2 = (P20' -P)' Ep = ^p2 + m2, 

(3.10) 
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and 4>s,d{po, |p|) are the partial scalar amplitudes, related to the corresponding partial 
vertices as 



0s,d(po, IpI) 



GsMpo^ IpI) 



M, 



D 



(3.11) 



-^P -Pt> 



Mr] is the deuteron mass, and the normalization factor reads M 



The components of the polarization vector of a vector particle moving with 4-momentum 
p = {E, p), having the polarization projection Ai = ±1, and mass M are 



e 



M 



Pi 



M 



,^M+P- 



Pi 



M 



where ^m is the polarization vector for the particle at rest. 



(3.12) 



^^' - "71 



/i\ 

i 



^--T. 



( I \ 



V " ; 







(3.13) 



The Dirac spinors, normalized as u{p)u{p) = 2m and v{p)v{p) = —2m, read 

f . \ 



u[p,s) 



m + Er, 



Xs 

crp 



Xs 



v[p,s) 



m + Er 



( "P ^^\ 

m + Ep 



\m + Ep^ ) 



V 



Xs 



(3.14) 



/ 



where Xs = ~'i-o'yXs, and Xs denotes the usual two-dimensional Pauli spinor. In gen- 
eral, the BS amplitude consists on eight partial components. As already mentioned, in 
eq. (3.9) we take into account only the most important ones, namely the S and D par- 
tial amplitudes. The other six amplitudes may become important at high transferred 
momenta [37, 38], hence for the present near-threshold process (3.1) they may be safely 
disregarded. Substituting eqs. (3.9)-(3.14) into (3.8) one obtains after some algebra [21] 



—t 



\M, 



1 



V 07I 

rfV Gs - Gd/V2 



3i 



„,^ 2Ep(27r)4(Mz,/2-Epj--p5 



A 



Ml 



+ 



E 



d'^p 



2 _ 2^^Sp,s„;rr 

Gd 



VTO^j^J 2Ep{2TiY {Mn 12 -Epf-pl 



'(Pl, P2; Pi, P2, ^v)5r+r'-Mcl (3-15) 

^t'Jn;rr'(Pi> P2; v'iV'2^4>)xt' {(^'ri)Xr {n^M^). 
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By closing the integration contour in the upper hemisphere and picking up the residuum 
in pq = Md/2 — Ep and introducing the notion of the deuteron S and D wave functions 
as 

. X ^5(^0, |p|)/47r GdIpo, |p|)/4vr 



v/2M^(2Ep - Md) ' '"^^ V^ME{2Ep - M^) ' 

2 r 
with - \p\'^d\p\{ul + ul)^l, (3.16) 

TT J 

the final expression for the amplitude T may be cast in the form 






An ^J E,{27i] 



\M, 



us{P) 



Sr+r';Md-^-^7K-{^*Md'^) Xr'(crn)Xrk (3.17) 



V2 

where n is a unit vector parallel to p. Equation (3.17) may be written in a more famil- 
iar form as to better emphasize the relation of our formulae with their non-relativistic 
analogues. For this sake observe that in eq. (3.17) one has J\M.d\ + ^^r+r'-Md — 



-v^(i rlr'llMd) and U*Mn) (crn) = (47r/3) E(-l)"+''e«'T-/3Yi«(n)Yi^(n), where 
Yiin(r7.) are the usual spherical harmonics. Then, by making use of the addition theo- 
rem for the spherical harmonics Yi^{n) and the Wigner-Eckart theorem for the matrix 
elements between states with definite angular momenta, the amplitude T^^'^-'^'^{s,t) be- 
comes 

Ks'^'^'is^t) = ^oY.j ;^4^<,L..'(Pi,P2;P,-P,Py) X (3.18) 



1 1 



'\lMd)Y,^{n)us{p)- J2 {2mliy,2\lMd)(^\r^r'\liy,2)Y;jn)unip)\ 



where the last line exactly coincides with a non-relativistic spin overlap between the 
deuteron wave function and two Pauli spinors of intermediate nucleons. Note that usu- 
ally in non-relativistic meson-nucleon theories the analogue of the amplitude T;^^-^''(s, t) 
is obtained by a non-relativistic reduction of the initially covariant operator O with sub- 
sequent calculations of the matrix element 



'SlS2 

2 



T^MvMd ^^■^■(^s,t) = {D,Md 0''-^-(12; 1'2V) 



\N.R. 



Pl,Sl,P2,S2 



Y, {D,Md \h,r,k2,r'){h,r,k2,r' |6^-^-(12; l'2V) |pi, Si, p2, ^2), (3.19) 



r,r'=l 
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where now (D, Aid \ki, r, k2, r') is the non-relativistic (complex conjugated) deuteron wave 
function in the momentum space. In our case, the use of the BS formahsm allows to com- 
pute the matrix element (3.18) directly from the covariant expression of the previously 
found operator O, avoiding the cumbersome procedure of its non-relativistic reduction for 
non-relativistic calculations by eq. (3.19). However, since in the considered kinematical 
range the BS wave functions, eq. (3.16), are rather similar to their non-relativistic ana- 
logues (see, e.g. [37]) and due to the formal equivalence of eqs. (3.18) and (3.19) we can 
use the former one for non-relativistic calculations as well by merely replacing us{p) and 
Uo{p) with the corresponding quantities from the Schrodinger equation with, e.g., Bonn 
or Paris potentials. 

B. Results for pn — > dV 

In our calculation of the cross section of the process (3.1) we use the numerical solution 
of the BS equation [25] in ladder approximation obtained with a realistic OBE interaction. 
The effective parameters used in the ladder approximation when solving the BS equation 
have been fixed in such a way to obtain a good description of the A^A^ elastic scattering 
data and the main static properties of the deuteron [37]. The obtained parameters turn 
out to be very close to those obtained in the non-relativistic framework of the Bonn group 
in determining the OBE A^A^ potential [26]. In this sense our analysis of the processes 
with the deuteron is consistent with the previous consideration of A^A^ -^ NNV processes. 

In figs. 11 and 12 the angular distribution and the total cross section are depicted. 
The shape of our evaluated cross section is rather similar to one computed in ref. [19]. 
However, there is some difference in the absolute values in case of 0. This difference 
is probably due to the fact that we fitted our parameters to vrA^ -^ (pN data and to 
A^A^ -^ (f)NN data, in particular, the new DISTO data point with a = 0.17 fib [6]. In 
addition, the methods of calculating the FSI effects are slightly different, which could 
provide some difference in values of FSI corrections as well. Nevertheless, the magnitude 
of our parameters and those of ref. [19] are very close to each other. It is also worth 
emphasizing that, as discussed in refs. [19, 28], there can be several sets of parameters 
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equally well describing the pp -^ ppV data. These sets differ not only by absolute values 
of parameters but also by the relative contributions of bremsstrahlung and conversion 
diagrams which possess different isospin structure and, hence differently contribute to pp 
and pn processes. So, in case of pn -^ dV processes the isospin transition corresponds 
to AT = 0, consequently the conversion diagrams are enhanced by a factor of "-3" 
in comparison with the bremsstrahlung diagram with exchange of neutral mesons. That 
means that (i) the contributions of bremsstrahlung terms are strongly reduced in pn -^ dV 
processes in comparison with the elementary pp reaction and (ii) the shape of the cross 
sections and angular distributions in the process pn -^ dV follows the behavior of the 
corresponding distribution in the elementary processes, modified by the deuteron wave 
function. This is clearly seen in fig. 11, where the shape of the angular distribution is 
very similar to the distributions found in pn reactions (the shapes of the corresponding 
distributions in figs. 2 and 5 must be compared with the curves labelled by open circles in 
fig. 11). At the threshold the distribution is fairly fiat, while with increasing energy some 
structure around the forward-backward directions becomes pronounced. The shape of 
the differential cross section in near forward and backward directions depends essentially 
upon the parameter set used in calculations (for a detailed analysis of the dependence of 
the shape of the angular distributions upon the chosen parameters, see [20, 22]). 

In fig. 12 the total cross sections for uo and meson production are depicted as a 
function of the excess energy. The uo production is larger than the cross section by 
roughly two orders of magnitude. From figs. 11 and 12 we observe that the computed cross 
sections exhibit a dependence on the potential used in the calculations of the deuteron 
wave function. This dependence is more pronounced at relatively low energies and has 
opposite behavior in uj (the BS wave function provides slightly larger cross sections than 
the Bonn one) and (j) production (the BS cross section is smaller than the Bonn one). 
This difference decreases with increasing energy. This can be observed from a comparison 
of angular distributions at two different energies (compare the curves labelled with open 
circles with the ones labelled by triangles in fig. 11). In fig. 12 this effect is seen for oo 
production (the upper panel) and not yet visible for mesons (lower panel). 

Eventually, in fig. 13 we analyze the values of the OZI ratio defined at equal excess 



20 



energies (dashed line) and at equal beam energies (solid line). Corresponding to our 
approach in both cases the FSI effects are the same (i.e., the deuteron in the final state). 
The only difference should appear due to the difference in the phase space volumes for 
uo and mesons which is estimated to be minimized for the ratio at equal excess energy 
near the threshold and more pronounced for the ratio at equal beam energies. As the 
energy increases the difference in the two definitions should not be too large, as seen in 
fig. 13. In both cases the OZI ratio is essentially larger than the one expected form the 
naive OZI rule restrictions. 

IV. SUMMARY 

In summary we present a combined analysis of vector meson production in pp -^ ppV, 
pn — * pnV and pn -^ dV processes for V = uj,(f). The elementary reactions A^A^ are 
treated within an effective meson-nucleon theory with most parameters fixed from inde- 
pendent experiments. The few remaining parameters are fitted to achieve a reasonably 
good description of the available angular distributions. With these parameters at our 
disposal we obtain a fairly good description of the energy dependence of total cross sec- 
tions. Then the OZI rule is analyzed for different definitions of the relevant ratio of (p/uj 
meson production. It is found that an enhancement of the OZI rule ratios can be ob- 
tained without any expected violation of the original rule and could be interpreted as 
dynamical effect which occurs as a sophisticated interference of different types of dia- 
grams and isospin effects. Using the same meson-nucleon theory we investigate vector 
production with a deuteron in the final state. In order to be able to use directly the 
results from elementary reaction we apply the Bethe-Salpeter formalism for the deuteron 
in the final state. The final expressions are presented in a form fully corresponding to 
the non-relativistic approach. Calculations with the same set of parameters show that 
within the proposed approach one can obtain reasonable values of the total cross section, 
close to the preliminary experimental data from ANKE. The OZI rule ratio is found to 
be almost independent of the energy and, similarly to the A^A^ case, is enhanced relative 
to the naive expectations, based on the OZI rule restrictions. 
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The proposed approach allows to calculate a number of polarization observables in the 
considered processes which could be directly related with the OZI rule. However this task 
is beyond of the goal of the present paper. An analysis of polarization observables in 
NN — > NN(f) reactions can be found in ref. [28] and a for pn — ^ d(f) in ref. [21]. 

Our approach extends the study of [22] by including the neutron channels and the 
deuteron final state. First experimental data for pn -^ du is at our disposal, while for 
pn -^ d(f) data are expected soon from COSY-ANKE. It should be noted that in [22] also 
a few nucleon resonances are taken into account for uj. This leads to a renormalization of 
parameters, coupling constants and cut-off's, and invokes further new parameters. The 
role of baryon resonances is, in particular, stressed in [42] . The comprehensive extension of 
our approach with inclusion of the full hst of resonances deserves a separate investigation, 
e.g. along the lines of [39-41]. 
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FIG. 1: Graphical representation of the scattering operator 0(12; \'2'V) (with truncated nucleon 
lines) for the vector meson production in A^A^ interactions within an effective nieson-nucleon 
theory. The dashed and solid thick lines symbolize the OBE of vr, o", w, /? mesons with effective 
parameters listed in Table I (see also ref. [26]), the wavy lines represent the produced vector 
meson. The first group of diagrams (a) corresponds to meson production from bremsstrahlung 
processes, while the last two diagrams (b) depict meson production from the internal conversion 
of the virtual up into real oj oi (j) mesons. The dots symbolize the dressing of the corresponding 
lines of virtual particles with monopole cut-off form factors. 
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FIG. 2: Angular distributions of uj mesons for the reaction NN -^ NNuj at the excess energy 
As 2 = 173 MeV. Dashed Unes correspond to contributions of conversion currents (fig. lb), 
sohd Unes represent the total contributions of bremsstrahlung and conversion diagrams (figs, la 
and lb). FSI is included in all contributions. Experimental data are from the COSY-TOF 
Collaboration [31]. 
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FIG. 3: Total cross sections of near-threshold production of to mesons in the process A^A^ — > 
NNuj (left panel: pp reaction, right panel: pn reaction) as a function of the energy excess Asa = 
■y/i — 2?7ip — m^j . Dashed lines correspond to contributions of conversion currents solely (fig. lb), 
dotted and solid lines represent results with including bremsstrahlung and conversion diagrams 
(figs, la and lb) without and with FSI effects taken into account, respectively. Experimental 
data are from SATURNE [32], COSY-TOF Collaboration [31] and DISTO Cohaboration [5, 6] 
(cf. the compilation in [22]). 
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FIG. 4: Ratios of the total cross sections of uj meson production in p-p and pn channels as a 
function of the excess energy. Dashed line denotes the ratio of cross sections without FSI; solid 
line reflects the ratio of total cross sections with FSI taken into account. 
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FIG. 5: Angular distributions of <^ mesons in the process NN -^ NN(j) (left panel: pp reaction, 
right panel: pn reaction) at the energy excess Asa = 83 MeV. Experimental data are from the 
DISTO Collaboration [5, 6]. Notation as in fig. 2. 
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FIG. 6: Total cross sections of near-threshold production of (j) mesons in the process A^A^ — > NNcj) 
as a function of the energy excess Asa = ^ — 2mp — m^. Experimental data is from the DISTO 
Collaboration [5, 6]. Notation as in fig. 3. 
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FIG. 7: Ratios of the total cross sections of (p meson production in pp and pn channels. Notation 
as in fig. 4. 



29 



o 




50 100 150 200 250 300 



AS^^^ [MeV] 



FIG. 8: Ratio of the total 0-to-u; production cross sections in pp and pn reactions. Dashed 

Unes denote the ratio of "net" cross sections, i.e. without any FSI; sohd hues reflect the ratio of 

total cross sections with FSI taken into account. The shaded area indicates the range of ratios 

I fiW j ^ 0.0071 and ( ^M^ j ~ 0.0079 (cf. Table I) used as input for the calculations. 
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FIG. 9: Ratio of the total 0-to-u; production cross sections in pp (solid lines) and pn (dashed 
lines) reactions as a function of the excess energy above the (p threshold. Left (right) panel 
corresponds to ratio of total cross sections with (without) FSI effects taken into account. Data 
for the PP reaction is from DISTO Collaboration [6]. 
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FIG. 10: Diagrams contributing to the process pn — > dV. The sum of bremsstrahlung (a) and 
conversion- type (b) diagrams results in the matrix element of the operator 0(12; 1'2'y) (fig. 1) 
sandwiched between the initial two-nucleon states and the final deuteron. 
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FIG. 11: Angular distributions for lo (left panel) and (right panel) mesons produced in pn — > 
dV processes for various values of the excess energy As"*^' ^ = ^fs — M^ — my {V = uj, (j)). Solid 
lines: results of calculations with the deuteron wave function obtained within the BS formalism; 
dashed lines: the non-relativistic deuteron wave function with the Bonn potential. The curves 
labelled with open circles correspond to values of As^''^ for which the experimental data of 
elementary process pp exits (cf. figs. 2-6); the curves labelled with triangles are presented for an 
illustration of the change of the distributions with the excess energy. 
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FIG. 12: Total cross sections for uj and (j) mesons for pn -^ dV reactions as a function the 
excess energy As2 = y^ — Md — my {V = uj,(j)). Solid lines: results of calculations with the 
deuteron wave function obtained within the BS formalism; dashed (dotted) lines: using the 
non-relativistic deuteron wave function with the Bonn (Paris) potential. The experimental data 
are from COSY-ANKE [35]. 
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FIG. 13: The OZI rule ratio for pn -^ dV processes as a function of the excess energy. The 
dashed Hne represents the ratio at equal values of the excess energy, while the solid line reflects 
the ratio at equal beam energies and is measured from the threshold of (p meson production. 
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